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Abstract 

Contributions to QCD non-singlet forward evolution kernels P(z) for the DGLAP equation and 
to V(x, y) for non-forward (ER-BL) evolution equation are calculated for a certain class of diagrams 
which include renormalon chains. Closed expressions are presented for the sums of these contribu- 
tions that dominate for a large number of flavors Nf 3> 1. Calculations are performed in covariant 
£-gauge, in the MS scheme. The assumption of "naive nonabelianization" approximation for kernel 
calculations is discussed. The partial solution to the ER-BL evolution equation is obtained. 
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1 Introduction 



Evolution kernels are the main ingredients of the well-known evolution equations for parton distri- 
bution of DIS processes and for parton wave functions in hard exclusive reactions. These equations 
describe the dependence of parton distribution functions and parton wave functions on the renormal- 
ization parameter n 2 . Here, I continue to discuss the diagrammatic analysis and multiloop calculation 
of the forward DGLAP evolution kernel P(z) Q and non-forward Efremov-Radyushkin-Brodsky- 
Lepage (ER-BL) kernel V(x, y) p in a class of "all-order" approximation of the perturbative QCD 
that has been started in ||. There, the regular method of calculation and resummation of certain 
classes of diagrams for these kernels has been suggested. These diagrams include the chains of one-loop 
self-energy parts (renormalon chains) into the one-loop diagrams (see Fig. 1). In this letter, the results 

for both the kinds of kernels, obtained earlier in the framework of a scalar model in six dimensions 

Nf 

with the Lagrangian Li nt = ^ (V^Vw)^) with the scalar "quark" flavours (tpi) and "gluon" (</?), are 

i 

extended to the non-singlet QCD kernels. For the readers convenience some important results of the 
previous paper Q would be reminded. 

The insertion of the chain into "gluon" line (chain- 1) of the diagram in Fig.l a,b and resummation 
over all bubbles lead to the transformation of the one-loop kernel Pq(z) = az = a(l — z) into the kernel 



Po(z) 



az 



chain! 



P^{z-A) 



az 



(z)- A (l-A) 



7^(0) 



where A = aNfy^O), a 



(An) 



(1) 



Here, 7^,( ¥ ,)(e) are the one-loop coefficients of the anomalous dimensions of quark (gluon at Nf = 1) 
fields in D-dimension (D = 6 — 2s) discussed in ||; for the scalar model 7^(e) = 7^(e) = B(2 — e,2 — 
e)C(e), and C(e) is a scheme-dependent factor corresponding to a certain choice of an MS -like scheme. 
The argument A of the function 7„ (A) in ([!]) is the standard anomalous dimension (AD) of a gluon 
field. So, one can conclude that the "all-order" result in ([!]) is completely determined by the single 
quark bubble diagram. The resummation of this "chain- F' subseries into an analytic function in A 
shouldn't be taken by surprise. Really, the considered problem can be connected with the calculation 
of large Nf asymptotics of the AD's in order of 1/Nf. An approach was suggested by Vasil'ev and 
collaborators at the beginning of 80' Q to calculate the renormalization-group functions in this limit, 
they used the conformal properties of the theory at the critical point g = g c corresponding to the 
non-trivial zero g c of the D-dimensional /3-function. This approach has been extended by J. Gracey 
for calculation of the AD's of the composite operators of DIS in QCD in any order n of PT, |5|, ||]. 
I have used another approach, which is close to |?], ||; contrary to the large Nf asymptotic method 
it does not appeal to the value of parameters NfTn, Ca/2 or Cp, associated in QCD with different 
kinds of loops. Following this way, the "improved" QCD kernel P^(z;A) has been obtained in || for 
the case of quark or gluon bubble chain insertions in the Feynman gauge. 

In this paper, we present the QCD results similar to Eq.(|l]), for each type of diagrams appearing 
in the covariant £- gauge for the DGLAP non-singlet kernel P(z;A). The analytic properties of the 
function P(z; A) in variable A are analyzed. The assumption of "Naive Nonabelianization" (NNA) 
approximation || for the kernel calculation |1C] is discussed and its deficiency is demonstrated. The 
ER-BL evolution kernel V(x, y) is obtained in the same approximation as the DGLAP kernel, by using 
the exact relations between P and V kernels [11, ]3| for a class of "triangular diagrams" in Fig. 1. The 
considered class of diagrams represents the leading Nf contributions to both kinds of kernels. At the 
end, a partial solution for the ER-BL equation is presented (compare with |L0], p^|). 

The obtained results are certainly useful for an independent check of complicated computer calcu- 
lations in higher orders of perturbation theory (PT), similar to |fl5|| ; they may be a starting point for 
further approximation procedures. 
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2 Triangular diagrams for the DGLAP evolution kernel in QCD 



Here, the results of the bubble chain resummation for QCD diagrams in Fig.l a,b,c for the DGLAP 
kernel are discussed. These diagrams generate contributions ~ a s (a s ln[l/z]) n in any order n of PT. 



+MC 






(a) 



(b) 



(c) 



+ 



+ 



Figure 1: The diagrams in figs, la - lc are the "triangular" diagrams for the QCD DGLAP kernel; 
dashed line for gluons, solid line for quarks; black circle denotes the sum of all kinds of the one-loop 
insertions (dashed circle), both quark and gluon (ghost) or mixed chains; MC denotes the mirror- 
conjugate diagram. 

Based on the resummation method of Ref. j| in the QCD version, one can derive the kernels 
p(ia,b,c) corresponding to the diagrams in Fig.l in the covariant £— gauge 



p( la \z;A) 
P (lb \z;A) 
P {lc \z;A) 



2z 



ci s Cf2z 
a s C F 2 ■ 
ci s Cf ■ 8(1 — z 



1 - Afz 



7,(0) 



l-A 



1 



7 9 (0) \ 



ci s Cf ■ 5(1 



1 



7 S (0) 



(l-A) lg (A) 



( A(3- 2A) 7fl (0) 



(2-A){l-A) lg {A) 



(3) 
(4) 



where a 9 



n s 



(47T) 



C F = (N 2 C - 1)/2N C , C A 
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N c and Tr = — are the Casimirs of SU(A r c ) group, 



and A = —a s j g (0). The function r ) g (e) is the one-loop coefficient of the anomalous dimension of gluon 
field in D-dimension, here D = 4 — 2e. In other words, it is the coefficient Z\{e) of a simple pole in 
the expansion of the gluon field renormalization constant Z, that includes both a finite part and all 



the powers of the e-expansion. Equations (|2|) - (Q) are valid for any kind of insertions, i.e., j g 
for the quark loop, 7 9 = 'jg 9 ^ for the gluon (ghost) loop, or for their sum 

lg (A,o=i g q) (A)+i { g 9) (A,0; 



7s 



(«) 



when both kinds of insertions are taken into account. The sum of contributions (J2j) , (||), (||) results in 
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P^ l '(z;A,^) which has the expected "plus form" 



a s P (z) 



a s C F 2 
(i s Cp1 



zz 



l (l-^) 2 + 



2z 



l-A 



7 g (0,0 



+ 



2z 



1 



(5) 
(6) 



where, for comparison, the one-loop result a s Po(z) is written down, the latter can be obtained as the 
limit pV)(z;A -» 0,0- Note that in © the 6(1 — z) - terms are exactly accumulated in the form 
of the [...]+ prescription, and the £ - terms successfully cancel. This is due to the evident current 
conservation for the case of quark bubble insertions, including the gluon bubbles into consideration 
merely modifies the effective AD, j g (A,£) — > j g 9 \A), conserving the structure of the result (||), see 
[||. Substituting the well-known expressions for j g (e) from the quark or gluon (ghost) loops (see, e.g., 

m 



7^M) 



-8N f T R B(D/2,D/2)C(e), 
^/2-^/2-l)((f^) + 

(l-^D-3) + (^) 2 e) C(e), 



(7) 



(8) 



into the general formulae @ - @ , and (|5|) one can obtain P^ (z; A, £) for both the quark and the gluon 
loop insertions simultaneously. Here, the coefficient C(e) = T(l — e)T(l + e) implies a certain choice 
of the MS scheme where every loop integral is multiplied by the scheme factor T(D/2 — l)(// 2 /4-7r) £ . 
The renormalization scheme dependence of P^ l \z; A) is accumulated by the factor C(e) 0. Of course, 
the final result (Jsj) will be gauge-dependent in virtue of the evident gauge dependence of the gluon 
loop contribution 7p (£,£), in this case, e.g., 



MO 



.7*(o,e) 



(0) + 7^(0,0 



(9) 



is the contribution to the one-loop renormalization of the gluon field. The positions of zeros of J g (A, £) 

in A, i. e., the poles of P(z; A, £), also depend on £. The kernel P^ (z; A) became gauge-invariant in the 

4. 



case when only the quark insertions are involved, i.e., 7 9 



19 ! 



A = A^ = -aA^(0) 



a s ^T R N f , 



and pW(z;Ato) -» P^\z;A) as it was presented in It is instructive to consider this case in detail. 
To this end, let us choose the common factor 7^ (0)/7g 9 ^ (A) in formula @ for the crude measure of 
modification of the kernel in comparison with the one-loop result a s Po(z). Considering the curve of 
this factor in the argument A in Fig. 2, one may conclude: 

(i) the range of convergence of PT series corresponds to the left zero of the function 7^ (A) and 
is equal to A = 5/2, that corresponds to a s o = Ihir/Nf, so, this range looks very broad f], a s < 5ir 
at N f = 3; 

(ii) in spite of a wide range of PT fidelity, the resummation into Pq 1 \z;A) is substantial - two 
zeros of the function Pq l \z;A) in A appear within the range of convergence (it depends on a certain 
MS scheme); 



2 For another popular definition of a minimal scheme, when a scheme factor is chosen as exp(c e), c = — 7b + . . . instead 
of r(D/2— 1), the coefficient C(e) does not contain any scheme "traces" in final expressions for the renormalization-group 
functions. 

3 Here we consider the evolution kernel P(z,A) (V(A)) by itself. We take out of the scope that the factorization scale 
fi 2 of hard processes would be chosen large enough, /1 2 > m 2 , where the p-meson mass represents the characteristic 
hadronic scale. Following this reason, the used coupling a s (/j, 2 ) could not be too large. 
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Figure 2: The curve of the factor r y g (0)/ / y g (A), the arrow on the picture corresponds to the point 
A = 1/ir, the first singularity appears at A = Aq = 5/2. 



(iii) the factor 7^(0) /j g g \A) decays quickly with the growth of the argument A. Really, if we 
take the naive boundary of the standard PT applicability, a s = 1 (at Nf = 3, A™ = 1/(2%)), then 
this factor falls approximately to 0.7 (at Nf = 6, = 1/tt it falls to 0.5, see arrow in Fig. 2); 

thus, the resummation is numerically important in this range. 

Note at the end that Eq.(|5|) could not provide valid asymptotic behavior of the kernels for z — > 0. 
A similar z-behavior is determined by the double-logarithmic corrections which are most singular at 
zero, like a s (a s ln 2 [z]J [14]. These contributions appear due to renormalization of the composite 
operator in the diagrams by ladder graphs, etc. rather than by the triangular ones. 



3 Analysis of the NNA assumption for kernel calculations 

The expansion of Pq 1 \z;A) in A provides the leading a s (a s Nf ln[l/z]) n dependence of the kernels 
with a large number Nf in any order n of PT Q. But these contributions do not numerically dominate 
for real numbers of flavours Nf = 4, 5, 6. That may be verified by comparing the total numerical 
results for the 2- and 3-loop AD's of composite operators (ADCO) in [15] with their Ay-leading terms 
(see ADCO in Table 1). Therefore, to obtain a satisfactory agreement at least with the second order 
results, one should take into account the contribution from subleading ./Vy-terms. As a first step, let 
us consider the contribution from the completed renormalization of the gluon line - it should generate 
a part of subleading terms. Below we shall examine two special choices of the gauge parameter £. To 
facilitate the diagrammatic analysis, it is instructive to inspect first the Landau gauge £ = 0. Indeed, 
the self-energy one-loop insertions into the quark lines as well as a certain part of vertex corrections to 
triangular diagrams are proportional to £; therefore, they disappear in the Landau gauge. Moreover, 
one should not consider the renormalization of parameter £. The analytic properties of the function 
P^(z;A,£ = 0) in the variable A = A(0) are modified - the function has no singularities in A until 
the "asymptotic freedom" exists, i.e., A < (at 13C.4 > ANf). In spite of all these profits the kernel 
(z; A, 0) generates the partial kernels a^Pn\(z), a^P^)(z), . . . which are rather far from the real 
ones. The ADCO 7(i,2) ip) corresponding to these kernels (here j(n) = Jq 1 dzz n P(z) ) are presented 
in Table 1. 
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Table 1. 



The results of 7(^2) ( n ) calculations performed in different ways, exact numerical results from [15], 
approximations obtained from P(z, A, £) with £ = and £ = —3; both numerical and analytical exact 
results are emphasized by the bold print. 





7(1) (") 


7(2) (n) 


CfCa 


Nj • Cp 


C\Cp 


N f ■ C F C A 


Nj ■ C F 


n=2 
Exact 

£ = -3 
£ = 


13.9 

11.3 
7.6 


64 
~27 


86.1 + 21.3 C(3) 

-42.0 
-13.2 


-12.9 - 21.3 C(3) 

12.9 
7.5 


224 
243 


n=4 
Exact 

p — _3 
£ = 


23.9 

23.5 
15.8 


13271 


140.0 + 19.2 C(3) 

-76.0 
-23.5 


-18.1 - 41.9 C(3) 

23. 
12.4 


384277 


2700 


243000 


n=6 
Exact 

£ = -3 
£ = 


29.7 

31.1 
20.7 


428119 


173 + 19.01 C(3) 

-95.6 
-29 


-20.4 - 54.0 C(3) 

28.5 
15.2 


80347571 


66150 


41674500 


n=8 
Exact 

£ = -3 


33.9 

36.3 
24 


36241943 


196.9 + 18.98 C(3) 
-109.0 
-33.0 


—21.9 — 62.7 C(3) 

32.3 
17.2 


-2.1619 


4762800 


n=10 
Exact 

£ = -3 

£ = 


37.27 

41.00 
27.29 


-8.5095 


216.0 + 18.96 C(3) 

-119.28 
-36.0 


-23.2 - 69.6 C(3) 

35.24 
18.68 


-2.3366 



Another exceptional gauge is £ = —3. For this gauge the coefficient of one-loop gluon AD 7 5 (0, —3) 
coincides with the coefficient bo of the /3-function 0. Therefore this gauge may be used for a refor- 
mulation of the so-called || NNA proposition to kernel calculations. To obtain the NNA result in 

a usual way, one should substitute the coefficient &o for 7^(0) into the expression for by hand 
(see, e.g., Jl(|). Note, the use of such an NNA procedure to improve Pq\z;A) leads to poor results 
even for a 2 s P\{z) term of the expansion; a similar observation was also done in [16|] . The NNA trick 

4 Here, for the /3(a s )-function we adapt (3(a s ) = — &o a s + • • ., bo = -5-CU — \^)'Tr 
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expresses common hope that the main logarithmic contribution may follow from the renormalization 
of the coupling constant. This renormalization appears as a sum of contributions from all the sources 
of renormalization of a s , corresponding diagrammatic analysis for two-loop kernels is presented in 
[0, |J|. In the case of the £ = — 3 gauge the one- loop gluon renormalization "imitates" the contri- 
butions from these other sources and the coefficient &o appears naturally. The elements of expansion 
of the ADCO j(n;A, — 3) (that corresponds to P'^(z;A, -3)) in a power series in a s , a? s 7(1) (n); 
a s 7(2) ( n )j ••• and a few numerical exact results from are collected in Table 1, let us compare 
them: 

(i) we consider there the contribution to the coefficient 7^) (n) which is generated by the gluon loops 
and associated with Casimirs CfCa/2, the C F -teim is missed, but its contribution is insignificant. 
It is seen that in this order the CpCA-terms are rather close to exact values (the accuracy is about 
10% for n > 2) and our approximation works rather well; 

(ii) in the next order the contributions to 7(2) (n) associated with the coefficients Nf ■ CfCa and 
C\Cf are generated, while the terms with the coefficients C F , Nf ■ C F , C f Ca are missed. In the 
third order, contrary to the previous item, all the generated terms are opposite in sign to the exact 
values, and the "£ = —3 approximation" doesn't work at all. So, we need the next step to improve 
the agreement - to obtain the subleading iVj-terms by the exact calculation. 

In any case, it seems rather difficult to collect the renormalization constant required by the NNA 
approximation in the kernel calculations. It is because different sources of renormalization of a s 
provide different z-dependent contributions, compare, e.g., Exp.© with Eq.(10) in g, the latter 
being generated by the insertions of self-energy quark parts into the quark line (chain 2). For this 
reason, necessary cancellation between the terms from different sources looks unlikely. 



4 Triangular diagrams for the non-forward ER-BL evolution kernel 



Here we present the results of the bubble resummation for the ER-BL kernel V(x,y). It can be 
obtained as a "byproduct" of the previous results for the kernel DGLAP P(z), i.e., in the same 
manner as it was done for the scalar model in Q. We shall use again the exact relations between the 
V and P kernels established in any order of PT [|llj for triangular diagrams. These relations were 
obtained by comparing counterterms for the same triangular diagrams considered in "forward" and 
"nonforward" kinematics. 

Let the diagram in Fig. la have a contribution to the DGLAP kernel in the form P(z) = p(z) + 
5(1 — z) ■ C; then its contribution to the ER-BL kernel is 



V(x,y)=C[9(y>x) 



-dz \+5(y-x)- C, 



(10) 



where C = 1 + (x — > x, y — > y). From relation ( |l0| ) and Eqs. (^), (|j) for p( lct ' c ) we immediately derive 
the expression for the sum of contributions l/( la+lc ) , 



V^ la+lc \x,y-A,0 =a s C F 2-C 



6(y>x)(l-A) 



x 



l-A 



1 



(l-A) 
(2 -A) 



7 (0,0 
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that may naturally be represented in the "plus form" . Expression (|ll| ) can be independently verified 
by other relations reducing any V to P ]TT| , [T7[ | (see formulae for the V — > P reduction there) and we 
came back to the same Eqs.@, @ for pT 1( w . Moreover, the first terms of the Taylor expansion of 
A) in A coincide with the results of the two- loop calculation in fLT| . The relation P — > V 
similar to Eq. (|To|) has also been derived for the diagram in Fig. lb 



V^ lb \x,y)=C 



(12) 
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therefore, substituting Eq.(||) into (12) we obtain 



V^ lb \x,y;A,0 = a s C F 2-C 



0(y >x)(~ 

y 



l-A 



1 



y-x 



7 g (0,0 



(13) 



Collecting the results in (|IT| ) and ( |l3|) we arrive at the final expression for in the "main bubbles" 
approximation 



vM(x,y;A,£) = a 8 C F 2-C 



6(y > x) 



l-A 



1-A + 



1 



y-x 



7 9 (0,0 



(14) 



which has a "plus form" again due to the vector current conservation. The contribution in 
( |j~4| ) should dominate for iV*/ ^ 1 in the kernel V. Besides, the function V^\x, y; A, £) possesses an 
important symmetry of its arguments x and y. Indeed, the function V(x,y; A,£) = V^'(x, y; A, £) ■ 
(yy) 1 ~ A is symmetric under the change x <-> y, V(x,y) = V(y, x). This symmetry allows us to obtain 
the eigenfunctions ip n (x) of the "reduced" evolution equation JH 



y; A)i> n {y; A)dy = 7 (n; A)4> n (x; A), 



(15) 



^n(y; A) ~ (yy) d ^ A) -^ Ct [A \y - y), here d^A) = (D A - l)/2, D A = 4- 2A, (16) 

and djf,(A) is the effective dimension of the quark field when the AD A is taken into account; C^\z) 
are the Gegenbauer polynomials of an order of a. The partial solutions <&(x;a s ,l) of the original 
ER-BL-equation ( where I = ln(/z 2 //Xo)) 

(fi 2 d^+l3(a s )d aB )$(x;a s ,l)= f V« (x, y; A) $(y;a s ,l)dy 



(17) 

are proportional to these eigenfunctions ip n (x;A) for the special case (3{a s ) = 0, see, e.g. ||. 

In the general case P(a s ) ^ let us start with an ansatz for the partial solution of Eq,fli~7|), 
$ n (x; a s , I) ~ Xn(o-s, I) • "^n{x; A), and the boundary condition is Xn(a s , 0) = 1; $ n (x; a s , 0) ~ ip n {x; A). 
For this ansatz Eq.([l7|) reduces to 



(/i 2 ^ 2 +/3(a s )5 as )ln($ (x; a s ,l)) = 7(n; ^4). 



(18) 



In the case n = the AD of the vector current 7(0; ^4) = 0, and the solution of the homogeneous 
equation in (fL8l) provides the "asymptotic wave function" 



$ofo a s , I) = A) ~ ((1 - x)xf 1 ~ A \ 



(19) 



where A = — a s (/i 2 )7(0, £) and a s (fj?) is the running coupling corresponding to (3{a s ). A similar 
solution has been discussed in fTO] in the framework of the standard NNA approximation. Solving 



simultaneously Eq. (18) and the renormalization-group equation for the coupling constant a s we arrive 
at the partial solution <fr n (x; a s , I) in the form 



$ n (x, a s ) ~ Xnip 2 ) ■ i>n{x] A); where Xn(^ 2 



exp 



a s (M 2 ) j(n,A) 



da 



(20) 



Recently, a form of the solution ~ t[) n (x;A) with A = —a s bo has been confirmed in |l^] by the 
consideration of conformal constraints |ll| on the meson wave functions in the limit Nf S> 1. 
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5 Conclusion 



In this paper, I present closed expressions in the "all order" approximation for the DGLAP kernel 
P(z) and ER-BL kernel V(x, y) appearing as a result of the resummation of a certain class of QCD 
diagrams with the renormalon chain insertions. The contributions from these diagrams, P^Xz-A) and 
V^ 1 ) (z; A) , give the leading Nf dependence of the kernels for a large number of flavours Nf ^> 1. These 

"improved" kernels are generating functions to obtain contributions to partial kernels like ai n+1 ^P( n ) (z) 
in any order n of perturbation expansion. Here A ~ a s is a new expansion parameter that coincides 
(in magnitude) with the anomalous dimension of the gluon field. On the other hand, the method 
of calculation suggested in [|| does not depend on the nature of self-energy insertions and does not 
appeal to the value of the parameters NjTr, Ca/2 or Cp associated with different loops. This allows 
us to obtain contributions from chains with different kinds of self-energy insertions, both quark and 
gluon (ghost) loops. The prize for this generalization is gauge dependence of the final results for 
PW(z;A) and V^\z;A) on the gauge parameter £. 

The result for the DGLAP non-singlet kernel P^\z; -A(£), £) is presented in (||) in the covariant 
£-gauge. The analytic properties of this kernel in the variable a s are discussed for quark bubble chains 
only, and in the general case for two values of the gauge parameter £ = 0; —3. The insufficiency of 
the NNA proposition for the kernel calculation is demonstrated by the evident calculation in the third 
order in a s (see Table 1). 

The contribution V^(x, y; A(£),£) to the non-forward ER-BL kernel (|l4]) is obtained for the same 



classes of diagrams as a "byproduct" of the previous technique [18|. A partial solution (PQ) to the 
ER-BL equation is derived. 
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